Orbital currents and charge density waves in a generalized Hubbard ladder 
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We study a generalized Hubbard model on the two-leg ladder at zero temperature, focusing on a 
parameter region with staggered flux (SF)/d-density wave (DDW) order. To guide our numerical 
calculations, we first investigate the location of a SF/DDW phase in the phase diagram of the half- 
filled weakly interacting ladder using a perturbative renormalization group (RG) and bosonization 
approach. For hole doping 8 away from half-filling, finite-size density-matrix renormalization-group 
(DMRG) calculations are used to study ladders with up to 200 rungs for intermediate-strength 
interactions. In the doped SF/DDW phase, the staggered rung current and the rung electron density 
both show periodic spatial oscillations, with characteristic wavelengths 2/8 and 1/5, respectively, 
corresponding to ordering wavevectors 2kp and Akp for the currents and densities, where 2kp = 
7r(l — 8). The density minima are located at the anti-phase domain walls of the staggered current. 
For sufficiently large dopings, SF/DDW order is suppressed. The rung density modulation also 
exists in neighboring phases where currents decay exponentially. We show that most of the DMRG 
results can be qualitatively understood from weak-coupling RG/bosonization arguments. However, 
while these arguments seem to suggest a crossover from non-decaying correlations to power-law 
decay at a length scale of order 1/8, the DMRG results are consistent with a true long-range order 
scenario for the currents and densities. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10 Pm, 71.30.+h 



I. INTRODUCTION 

The possibility of finding new phases of matter is a ma- 
jor motivation behind the study of materials of strongly 
correlated electrons. A phase with staggered orbital cur- 
rents which was first considered theoretically in 1968i 
and then rediscovered two decades latesS* 3 ^^ has lately 
been the subject of a revival of interest, mainly due to 
recent proposals^ that the pseudogap region 9 , in the 
phase diagram of the cuprate high-temperature super- 
conductors may be characterized by this kind of order, 
either long-ranged (i.e., a true broken-symmetry state) j£ 
or fluctuating^ This orbital current phase is known vari- 
ously as the orbital antiferromagnet, staggered flux (SF) 
or d-density waveitt (DDW) phase; in this paper we will 
refer to it as the SF/DDW phase. The long-range ordered 
version of this phase breaks the rotational and transla- 
tional symmetries of the underlying Hamiltonian, in addi- 
tion to time reversal symmetry The fundamental experi- 
mental signature of the long-range-ordered SF/DDW sce- 
nario with ordering wavevector (tt, it) is an elastic Bragg 
peak at that wavevector in neutron scattering&ii The 
results of some recent neutron scattering experiments on 
underdoped YBa2Cu306+:r have been argued to be con- 
sistent with this scenario .iLiSii 3 . A different circulating- 
current broken-symmetry state with a current pattern 
that does not break translational symmetry has also been 
proposed for the pseudogap phase in the cuprates— 

It is important to acquire an understanding of what 
kinds of microscopic models can give rise to a SF/DDW 
ground state. For models of interacting fcrmions in two 
spatial dimensions this issue has been addressed by many 



authora 2 ! 3 ! 4 ! 5 ! 6 : 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 ! 23 ! 24 ! 25 using a variety 
of methods. While the two-dimensional case is the most 
relevant one for the cuprates, the behavior of strongly 
correlated electrons in two dimensions is still a subject 
that is marred by great controversy. In a few special 
model cases the fermion sign problem is absent and thus 
quantum Monte Carlo simulations can be done reliably, 
leading among other things to interesting findings with 
regard to SF/DDW order— and also more exotic current- 
carrying states^ However, in the vast majority of cases, 
and certainly for the models and parameter values that 
are expected to be most relevant for real materials, the 
available analytical and numerical methods are only ap- 
proximate, and their reliability is difficult to gauge. 

In contrast, the fact that many powerful methods ex- 
ist for one spatial dimension has enabled much solid 
knowledge to be established about strongly correlated 
one-dimensional systems, and it is hoped that some of 
this may also be relevant for the behavior of corre- 
lated electrons in two dimensions. In this regard, two- 
leg Hubbard and t-J ladders have attracted much in- 
terest, as it has been found that these models have a 
spin gap and that upon doping away from half-filling 
the dominant correlations are d-wavc-likc superconduct- 
ing correlations^ 7 - two features which are reminiscent of 
the pseudogap and the d x 2_ y 2 symmetry of the super- 
conducting state in the cuprates, respectively. Ladder 
systems are also interesting in their own right, not least 
due to experimental realizations of such materials. The 
most well-known example, Sri4_ a; Ca a ;Cu2404i, contains 
two-leg ladder substructures, and has been found to have 
a spin gap, and to become superconducting upon doping 
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at high pressure i2L2& 

The possibility of having SF/DDW order in 
two-leg ladders has also received much attention 

r:c::n tl--. 29.30.31.32.33.34.35.36.37.38.39.40.41 rpj^ two-leg 

ladder is the simplest geometry that can support a 
SF/DDW phase, with currents flowing around the 
elementary square plaquettes. In this paper we present 
an extensive study of SF/DDW order in generalized 
Hubbard ladders, focusing mainly on the case of finite 
hole doping away from half-filling. Some of the results 
that will be discussed have been briefly presented in Ref. 

m 

We first use bosonization and a perturbative 
renormalization-group (RG) approach to identify a pa- 
rameter region with long-range SF/DDW order in the 
weakly interacting half-filled ladder i 3 ^ 7 ^ Finite-system 
density-matrix renormalization-group (DMRG) calcula- 
tions are then used to study the SF/DDW phase in doped 
ladders, for intermediate-strength interactions and lad- 
der sizes up to 200 rungs. For the (rational) hole dopings 
considered, currents are found to be large (of order the 
nearest-neighbor hopping amplitude t) and show no ev- 
idence of decay. As the doping is increased the currents 
in the SF/DDW phase decrease in magnitude, and for 
sufficiently strong doping SF/DDW order becomes com- 
pletely suppressed. The SF/DDW phase has a spin gap, 
but ci-wavc-like superconducting correlations decay expo- 
nentially. 

In the half- filled SF/DDW phase, the sign of the rung 
current changes from one rung to the next, while the 
magnitude stays the same, corresponding to an order- 
ing wavevector 2kp — 7r in the direction parallel to 
the legs. In contrast, the rung currents found in the 
doped SF/DDW phase vary both in sign and magni- 
tude, corresponding to an ordering wavevector 2kp = 
7r(l — 5) =/= 7r, where 5 is the hole doping away from 
half-filling. Such phases are often referred to as incom- 
mensurate. An incommensurate SF/DDW phase has 
been suggested as a candidate for the hidden order in 
the heavy fermion compound URu2Si2^ Incommensu- 
rate SF/DDW phases have also been studied in various 
two-dimensional modelsiiSiSM 3 . [A note on terminology: 
In the remainder of this paper only ordering patterns 
with an infinite unit cell will be characterized as 'incom- 
mensurate;' orderings with a finite unit cell (even when it 
is very large) will be referred to as 'commensurate.' The 
former (latter) case includes orderings with wavevector 
2kp = 7r(l — 5) with 5 an irrational (rational) number.] 

Charge density wave order in ladder systems 
is also a topic that has attracted much inter- 
est recently, both theoretically 37 i 38 i 44 i 45 i 46 i 47 i 48 and 
experimentallyi 49 i 50 i 51 i 52 i 53 i 54 Interestingly, we find that 
a charge density modulation coexists with the orbital 
currents in the doped SF/DDW phase. Like the cur- 
rents, this density modulation is not found to decay. 
The staggered rung current and the rung electron den- 
sity both show a periodic spatial variation, with charac- 
teristic wavelengths 2/5 and 1/5, respectively. The den- 



sity modulation corresponds to an ordering wavevector 
Akp, with two doped holes (one per leg) per wavelength. 
Furthermore, the minima of the electron density are lo- 
cated at the zeros of the staggered rung current. In other 
words, the rungs at which the doped holes predominantly 
sit are antiphase domain walls for the current pattern. 
This property, and the related factor of two ratio be- 
tween the periods of the staggered rung current and the 
rung electron density, are very reminiscent of the prop- 
erties of stripe phases with coexisting spin and charge 
order in doped antiferromagnets^ 5 " in which the role of 
the current in the SF/DDW phase is played by the spin 
density. The rung density modulation is found to per- 
sist in the neighboring phases of the SF/DDW phase in 
the phase diagram in which current correlations decay 
exponentially. We note that coexistence of current and 
charge order for an SF/DDW phase in two dimensions 
with ordering wavevector q ^ (jt, w) for the currents has 
previously been discussed within the framework of a phc- 
nomenological mean-field approach^ 3 " 

Is the current and charge density order long-ranged? 
The DMRG calculations do seem to support a commen- 
surate true long-range order scenario. Let us first note 
that, to the best of our knowledge, such a scenario does 
not violate any exact theorems that are directly applica- 
ble to the lattice model studied here; in particular, the 
Mcrmin- Wagner theorem is respected. Of course, these 
nonviolation arguments do not say whether or not true 
long-range order actually will occur. A weak-coupling 
RG/bosonization analysis of the continuum limit of the 
lattice model predicts that (1) true long-range order re- 
quires Umklapp interactions to be relevant in the RG 
sense, (2) at half-filling, Umklapp interactions are rel- 
evant, resulting in true long-range order j 3 ^ 7 ^ (3) for 
dopings away from half-filling Umklapp interactions are 
irrelevant, and as a result correlation functions decay as 
power laws (for distances much larger than 1/521), i.e., 
only quasi long-range order is predictedi 3 ^ 7 - DMRG cal- 
culations find true long-range order at half-filling^ thus 
agreeing with the RG/bosonization analysis in that case. 
But as already noted, DMRG calculations seem to sup- 
port a true long-range order scenario also for the doped 
ladder. We have not been able to resolve this appar- 
ent discrepancy between the predictions of the DMRG 
and RG/bosonization methods, and can only speculate 
about its origin. Both methods have their limitations 
and weaknesses. One particular question mark associ- 
ated with the weak-coupling RG/bosonization method 
is whether its predictions regarding the range of the or- 
der may be unreliable for the ladders studied here with 
DMRG, for which the interactions are not weak. But one 
certainly cannot exclude the possibility that there may 
be a different and perhaps quite subtle reason for the dis- 
crepancy. We refer the reader to the concluding Sec. IVII 
for an attempt at a fuller discussion of these issues. 

The paper is organized as follows. The generalized 
Hubbard model on the two-leg ladder is introduced in 
Sec. [H] In Sec. IIIII we use a perturbative RG ap- 
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proach and bosonization to study the phase diagram of 
the model for weak interactions at half-filling. In Sec. I1VI 
we present an extensive DMRG study of the model for 
nonzero doping and intermediate-strength interactions. 
In Sec. the doped ladder is discussed using weak- 
coupling RG and bosonization arguments. Sec. IVII con- 
tains a summary and concluding discussion. Some de- 
tails of the RG/bosonization calculations are discussed in 
three appendices. Throughout the paper the main focus 
is on the SF/DDW phase, but properties of neighboring 
phases in the phase diagram are also discussed. 



II. MODEL 

The model we will study is an "extended" or "general- 
ized" Hubbard model with various nearest-neighbor in- 
teractions in addition to the on-site Hubbard term. For 
the purpose of finding SF/DDW order, it is necessary 
to allow for these additional interactions because the 
ground states of the pure Hubbard ladder and the re- 
lated t-J ladder have been found to only have short- 
ranged SF/DDW correlationsi2ii2& Another, more gen- 
eral, motivation for studying generalized Hubbard mod- 
els is that, although much attention has been paid over 
the years to the two-dimensional Hubbard and t-J models 
as purportedly "minimal" effective models for describing 
the physics of the two-dimensional Cu02 planes in the 
cuprate superconductors, there is an increasing amount 
of studies that suggest that these models do not support 
high-temperature superconductivityi^Si 

The Hamiltonian of our model is H = Hq + Hj, where 
the kinetic energy and interaction operators are, respec- 
tively, 

#o = ( E c W c ^ + c t-C 6 x 8 +H.c, 

xs ^ I ' 

Hi = 2J ( [ Un £,x,1 n £,x,i + V\\nt x nt,x+i 

+ J\\Se x • Se yX+ i] + V±n tx nbx + J±S tx • St, x ^J . (2.1) 

The operator c\ xs creates an electron on rung x = 
1,...,L on the "top" or "bottom" leg £ — t,b with spin 
s =T,4, and obeys {q xs , c\, x>a) } = 6u'8 XX 'S ss >. Electrons 
can hop between nearest-neighbor sites along legs and 
rungs with hopping amplitude t. Electron density opera- 
tors arc given by n (xs = c\ xs ci xs and n gx = n £xs, and 

the spin operator is Sg x = jEss'is* 788 '^ 8 '' where 
the a 1 are Pauli matrices. The interactions consist of the 
Hubbard on-site term with strength U, as well as nearest- 
neighbor density-density and spin exchange interactions 
with strength V± and J± along the rungs, and Vj| and 
J|i along the legs. As the leg interactions will be seen to 
favor other phases than SF/DDW, we will in most of the 
paper take Vj| = J|| = 0. 



III. WEAKLY INTERACTING ELECTRONS ON 
THE HALF-FILLED LADDER 

In this section we study the case of weakly interacting 
electrons on a half-filled ladder i^ ^LS^ Th e fermionic 
perturbative RG approach of Refs. I58l59l60l is employed 
to obtain a low-energy effective theory which is then an- 
alyzed using bosonization and semiclassical considera- 
tions. Our focus is on finding a parameter regime for 
which the model (|2.1() has a ground state with long- 
ranged SF/DDW order. The phase diagrams obtained 
for the weakly interacting half-filled ladder serve as help- 
ful guides in the search for SF/DDW order in the doped 
ladder for intermediate interaction strengths, using the 
DMRG method, as described in Sec. IIVI Further- 
more, some of the results established here are used in the 
RG/bosonization analysis of the doped ladder in Sec. [V] 

A. Continuum description 

To study the weak-interaction problem, wc first diago- 
nalize the kinetic energy. As it is convenient to deal with 
a translationally invariant Hamiltonian, we impose pe- 
riodic boundary conditions in the direction parallel to 
the legs. By introducing even and odd combinations 
c\xs = {<hxs + (-1) A Qm;s)/\/2 with A = 1,2, and Fourier 
transforming along the leg direction, the kinetic energy 
is diagonalized in momentum space, describing an anti- 
bonding band (A = 1) and a bonding band (A = 2), with 
dispersions £\{k) = — 2t cosfc— (— l) x t (we set the lattice 
constant equal to 1). The Fermi level of the noninteract- 
ing half-filled system is at zero energy and crosses both 
bands, thus giving rise to four Fermi points ±/c.fi,2 which 
satisfy 

fcpi + k F2 = 7T. (3.1) 

As interactions are assumed to be weak, we focus on 
states very close to the Fermi energy, and linearize the 
kinetic energy around the Fermi points. The band oper- 
ator c\ xs can be expressed as a sum of right-moving (R) 
and left-moving (L) components, 

^^^e'"" 1 ^^), (3.2) 
p 

where P = R, L = ±1. The continuum fields ipp\ s (x) are 
slowly varying on the scale of the lattice constant. The 
linearized kinetic energy can be written Hq = J dxTlo , 
where 

H Q = -iv F E P iyp\ s dxil>p\s, (3.3) 

PXs 

where the Fermi velocity is vp — 

Interactions will scatter electrons between the Fermi 
points. Introducing the "currents' — 

J\p = ^i'pXs^PXs, J\p = 2 X]^As CTss '^ FAs '' 

s ss' 



4 



Lp = ^2 i'PU^Pis, L P = i ^2 i 3 Pl s a 'ss'i>P2 s ' , 
s ss' 

Mxp = -iippxfippxi, N Paa > = ippuipP2s', (3.4) 

the Hamiltonian density that describes the effects of the 
weak interactions to leading order can be written Hi = 
nf ] +Uf\ where 

T^P = —gip JxrJxl — 9xp(JirJ2l + J2rJ\l) 

A 

— 9lo ^ J\R ■ J\L — 9xo{JlR ■ J2L + J2R ■ Jil) 

X 

- 9t P {L R L L + L R L[) - g ta (L R ■ L L + L R ■ L[) 7 (3.5) 
Uf = -g xu {M lR Ml L + M 2R M\ l ) 

- Y.^iN Rss ,Nl ss , + g tu2 N Rss ,Nl s , s ) +11.C. (3.6) 



Here Ji \ represents (interband) Umklapp interactions. 



B. Bosonization 

Bosonization proves very helpful for interpreting the 
low-energy effective theory resulting from the perturba- 
tive RG flow (see Sec. IHI Cjl . In the Abclian bosonization 
formalismSiiSiSiSM^ the fermionic field operators ippxs 
can be expressed in terms of dual or conjugate Hermitian 
bosonic fields 4>\ s and 6\ s as 







f 



PXs 



V2i 



--Vx s exp[i(P<f> Xs + e Xs )}, (3.7) 



where a is a short-distance cutoff. The only equal-time 
nonzero commutator between the bosonic fields is (taking 
the limit aj{x — x') — > 0) 

[<l>x s (x),0x> s >(x')} = iTrdxx'S ss >H(x-x'), (3.8) 

where H(x) is the Heaviside step function. The long- 
wavelength normal-ordered fermionic densities can be ex- 
pressed in terms of the bosonic fields as 



V'pAsV'PAs = —d x {4> Xs + POxs)- 
Z7T 



For later use we also define 



NpXs = 



dx ippx s ippxs- 



(3.9) 



(3.10) 



The (Majorana) Klein factorsS*Si r]x s in Eq. (|3~7|) com- 
mute with the bosonic fields, and satisfy {ijxs, Vx's'} = 
ZSxx'Sss'- The Klein factor conventions used for the 
Hamiltonian and the order parameters considered in this 
paper (see below) are explained in Appendix 1X1 

Charge (p) and spin (a) operators are defined as 4>x P = 

(0AT + 4>xi)/V2, (px a = (0at - 0ai)/\/2. It is convenient 
to make a further change of basis, to <j) rv = (r4>\ v + 



02i^)/\/2, where r = ± and v = p,a. Similar definitions 
apply to the 0-operators. In the (r, z/)-basis, the kinetic 
energy density reads 



= if E [(^<M 2 + {d x e rv ) 2 ] . 



(3.11) 



Furthermore, H ( p = nf a) +H { j lb} , where 



TV 

W/ lb) = (9^)2 { 2 cos20 +CT g la cos 2<f>- a 
+ 9xa cos 26- a — g tlJ cos 2#_ p — cos 29- p 
{g ta ~ 4p t p) cos20_ CT + (g ta + Ag tp ) cos 26L CT |, (3.13) 

with g rp = -{gi p + rg xp ) and g r „ = -\{g\ a + rg xa ). 
Finally, the Umklapp interaction density reads 



(2) _ 



(2na)'- 



■ cos 2(f> 



g xu cos 29-p + gtui cos 2(f)- 



(ffiui + gtu2) cos2</) +(T + g tu2 cos26L CT . (3.14) 



C. Perturbative renormalization group approach 

The "bare" values of the nine coupling constants gi 
in the continuum field theory can be expressed in terms 
of the interaction parameters in Eq. I|2.1|l . and equa- 
tions describing the RG flow (to one loop) of these cou- 
pling constants near the noninteracting fixed point can 
be derived^ (see Appendix lB|l . These RG equations are 
then solved numerically for given sets of initial condi- 
tions. Typically, it is found that as the length scale is 
increased, some of the couplings grow (sometimes after a 
sign change), while other couplings remain small. 

The couplings which grow to be large (i.e., of order 1, 
at which point the numerical integration must be stopped 
since the perturbative RG equations are only valid for 
weak coupling) will tend to lock some of the bosonic 
fields in the minima of the cosine potentials in the Hamil- 
tonian. In fact, in the parameter region considered in 
this paper, at half-filling one of the conjugate fields 4> rl/ 
and 6 rv in each sector rv always becomes locked. These 
locked (and thus also gapped) operator fields can then 
be approximated by their c-number expectation values 
(integer multiples of 7r/2 for the Hamiltonian considered 
above), while their conjugate fields will fluctuate wildly. 
To deduce the nature of the ground state, it is necessary 
to bosonize various order parameters and use the infor- 
mation about the field-locking pattern to deduce which, 
if any, of these order parameters are nonzero, or have the 
most slowly decaying correlations. 
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cos 6-n sin#_ 



(3.17) 



SF/DDW 



— o — o — c> 

*T ^ s/ 

— o — o — o- 



FIG. 1: Schematic illustration of phases at half-filling in the 
strong-coupling limit. Open (filled) arrows represent currents 
(spins). The shading indicates singlet pairs. 



D. Phases, physical observables, and order 
parameters 

As our main interest here is the SF/DDW phase, we 
will focus on a parameter region where this kind of order 
is realized. Three other phases will also be seen to be 
present in this parameter region in the weakly interacting 
half-filled ladder. These are the CDW, D-Mott, and S- 
Mott phases & The nature of these four phases, most 
easily understood in the strong-coupling limit j28**i8<6Z<£& 
is illustrated schematically in Fig. ^ These phases are 
all completely gapped, i.e., both charge sectors ±p and 
both spin sectors ±er are gapped. Recently it has been 
shown that the two-leg ladder also can support four other 
completely gapped phases at half-filling^LS but none of 
these additional phases appear in the parameter region 
studied here. 

We now discuss in some detail the various phases and 
the order parameters/physical observables that charac- 
terize them, at half-filling^^i^^ (the bosonized form 
of these physical observables in the doped ladder will be 
discussed in Sec. IV A() . In the SF/DDW phase, cur- 
rents flow around the plaquettes as shown in Fig. ^ The 
ordering wavevector of the currents is (q x ,q y ) — (tt,tt), 
where q x (q y ) is the wavevector in the direction parallel 
(perpendicular) to the legs. The order parameter of this 
phase is 



OsF/DDW = COS(/) +p C0S(j) +r7 COS 0_ p COS#_ CT 

+ sin <f)+ p sin 4>+ a sin 6L P sin 6- a . (3.15) 

Denoting the rung current by j±(x), and defining the 
staggered rung current as 

j s (x) = (-l) x j±(x), (3.16) 

we have (j s (x)) cx (Osf/ddwOe)); note that the expecta- 
tion values are independent of x. 

In the CDW phase the electron density exhibits a 
checkerboard pattern. This charge density order is thus 
characterized by the same ordering wavevector (tt, tt) as 
the currents in the SF/DDW phase. The order parameter 



In the CDW phase the expectation value of the deviation 
6ni(x) of the charge density at leg £, rung x from the av- 
erage density is given by (Snt^ix)) cx ±(— l) x (Ocdw(x)}- 
Both the SF/DDW and CDW phase spontaneously 
break a Z2 symmetry, and thus have a two-fold degen- 
erate ground state with true long-range order. This is 
intimately related to the fact that the order parameters 
of these phases depend on <j) +p which is locked by the 
Umklapp interaction, Eq. I|3.14|l . In contrast, the D- 
Mott and S-Mott phases do not break any symmetry, 
but are characterized by a unique ground state with ex- 
ponentially decaying <i-wave-likc and s-wavc-like super- 
conducting (SC) correlations, respectively. The pairing 
operators can be taken to be 

AdSC/SScM = {^Rl^Lli+^Ll^Rli) 

T {lpR2\lpL2i+lpL2\^R2i)- (3.18) 

In Adsc the components of bands 1 and 2 (correspond- 
ing to transverse wavevector 7r and 0, respectively) come 
with opposite signs, while in Agsc they have the same 
sign; this justifies calling these operators d- and s- 
wave-like, respectively. One finds that (A a (x)At (0)) cx 
(O a (x)Ot(0)>, where a = DSC or SSC, and the order 
parameters are 



Odsc 
Ossc 



e ie+p (cos 4>+a COS 9-p COS <\)- a 

—i sin</> +cr sinfl-psini/i-o-), (3.19) 
, sin </)_ (T 

(3.20) 



e 46,+p (sin <j) +a cos i 



O 



CDW 



IJ +P ' 



sin tt-„ cos 



The exponential decay of the pairing correlations comes 
from the e* e+p factor, as 9 +p is conjugate to the locked 
field 4>+ p - 

The field-locking patterns that realize the SF/DDW 
and CDW phases are easily deduced from their respective 
order parameters. Due to the absence of (f>+ p in the DSC 
and SSC order parameters, the value of (4>+ p ) in the D- 
Mott and S-Mott phases must either be determined from 
the RG flow, or it can be deduced from the fact that 
the SF/DDW and the CDW phase can be regarded as 
the Ising-ordered counterparts of the quantum disordered 
D-Mott and S-Mott phase, respectively^ Field-locking 
patterns for the four phases are given in Table [I] 



E. Weak-coupling phase diagrams 

Fig. HJa) shows the phase diagram for Jii = Vji = 
and V± > in the weak-coupling limits For the region 
U > we may compare our results to those in Ref. l38l 
and find excellent agreement. The SF/DDW phase ap- 
pears between the D-Mott and CDW phases and is seen 
to extend well beyond the region where it was first re- 
alized to exist which was characterized by exact or 
approximate SO (5) symmetrySiS (the SO (5) symmetry 
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Phase 


(M 


(M 


(0- P ) 






SF/DDW 


0, TT 













CDW 


0, TV 





tt/2 







D-Mott 
















S-Mott 








tt/2 








TABLE I: Field-locking patterns (only unique up to a global 
gauge transformation of the fields) characterizing the ground 
states of the four phases encountered at half-filling. The two- 
fold degeneracy of the SF/DDW and CDW ground states has 
been reflected in the two values given for {<f>+ P ). As Q- a 
and (p-a are conjugate fields, it follows that if one of them 
is locked, the other is strongly fluctuating; this is indicated 
by 



occurs along the dotted line in the figure). The nature 
of the quantum phase transitions between the various 
phases has been discussed in detail in Refs. I3^l3§ (see 
also Ref. [EE 



Next, we consider the effects of including interactions 
along the legs. We fix U/V± = 0.25 and add a finite 
J\\/V± (Fig. 0») or VyVl (Fig. He)). J^/V± is seen 
to favor the D-Mott phase. This squeezes the SF/DDW 
phase which eventually disappears for large enough val- 
ues of J\\/V±. On the other hand, V\\/Vj_ is seen to favor 
the CDW phase, which is easy to understand intuitively. 
However, in this case the SF/DDW order is more ro- 
bust: the width of the SF/DDW region remains more or 
less constant as V\\/V± is increased, while its location is 
pushed towards larger values of J±/V±. 

From Fig. [2] the conditions that favor SF/DDW order 
for repulsive interactions and positive exchange constants 
can be roughly summarized as follows^ the ratios J±/U 
and V±/U must be sufficiently large, but not too large 
(as that favors the D-Mott and CDW phase respectively), 
with leg interactions Jm and Vji preferably zero or at least 
small compared to the other interaction parameters. One 
might speculate that the reason why anisotropic interac- 
tions are required to stabilize SF/DDW order in the two- 
leg ladder could be related to the fact that the currents 
in the SF/DDW phase must themselves be anisotropic 
in this geometry in order to satisfy current conservation 
(the rung currents are twice as large as the leg currents). 
It is also interesting to note that in mean-fielcUiSi and 
renormalization-group2L2£ studies of generalized Hub- 
bard models on the square lattice in two dimensions, at 
or near half- filling, and with isotropic interactions (i.e., 
J_i_ = J|| = J and V± = V\\ = V), the conditions that are 
most favorable for SF/DDW ordering tendencies also ap- 
pear to be that the ratios J/U and 4V/U be sufficiently 
large. 
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FIG. 2: Weak-coupling phase diagrams at half-filling, for 
V± > 0. Of the four phases (illustrated in Fig. that appear 
in these phase diagrams, the SF/DDW and CDW phase have 
true long-range order. The SF/DDW phase appears between 
the CDW and D-Mott phase, (a) Jy = Vj, = 0. The half-filled 
model has SO (5) symmetry along the dashed line. The phase 
boundary line U + (3/4) Jx = V± is independent of the filling 
and valid for arbitrary interaction strengths, as shown in Refs. 
liollill . (b) Vj| = 0, U/V± = 0.25. (c) J,| = 0, U/V± = 0.25. 
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FIG. 3: (Color online) In real-valued DMRG, possible long- 
ranged current patterns are detected by calculating current- 
current correlations centered around the ladder middle. 
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FIG. 4: (Color online) In complex-valued DMRG, possible 
long-ranged current patterns are detected by introducing a 
symmetry-breaking infinitesimal surface current and observ- 
ing whether it sustains a long-ranged current pattern. 



IV. DMRG STUDY OF THE DOPED LADDER 

In this section we discuss an extensive DMRG£2iZ2i£i 
study of the ground state of hole-doped two-leg lad- 
ders with up to 200 rungs for intermediate interaction 
strengths. The doping away from half-filling is measured 
by the parameter 5, defined as 

im !-» = §, (4.1, 

where n is the average number of electrons per site, and 
Mh and L are the number of (doped) holes and rungs, re- 
spectively. Open boundary conditions were used, as these 
are computationally preferable in DMRG. The weak- 
coupling phase diagrams at half-filling discussed in Sec. 
IIII El were used as guides to suggest interesting parameter 
regimes. 

To detect the existence of orbital currents, we have 
considered two approaches: 

(i) Within a standard DMRG calculation in real num- 
ber space, we have studied the decay of rung-rung cor- 
relations where the rungs have been centered about the 
middle of the ladder to minimize boundary effects (Fig. 
EJ): 

C(r) = (j± (L/2 + r/2)j± (L/2 - r/2)>, (4.2) 

where the rung current jj_(x) = J(t yX ),{b.x) an d the cur- 
rent operator between two nearest-neighbor sites i and j 
is given by 

''Xf'L.'j- 'j..-'-)- (4-3) 

(ii) Within a DMRG calculation generalized to com- 
plex number space, one may obtain non-vanishing cur- 
rent expectation values, 

Wl * 0. (4.4) 

In the textbook case of a ferromagnet that develops a fi- 
nite magnetization while spontaneously breaking the con- 
tinuous rotational symmetry, this is achieved in theoret- 
ical approaches by forcing a possible symmetry breaking 
by introducing an infinitesimally strong external fieldS 
In close analogy, we break the Z-i symmetry of the direc- 
tion of plaqucttc currents by introducing an infinitesimal 
surface current —hjj_(l) on the first rung of the ladder, 
corresponding to a weak surface field (Fig. 0J. To assure 



that this infinitesimal surface (or edge) current does not 
change the physics of the ladder under study, we have 
taken it down to values of O.OOOli (in some special cases 
even lower values were checked). We find that results 
do not change qualitatively while decreasing the surface 
current: in the phase that does support SF/DDW order, 
the magnitude of the orbital current measured numeri- 
cally docs not depend on its value in the small current 
limit, while in the phases that do not support SF/DDW 
order, the current decays exponentially fast into the bulk 
and its magnitude is proportional to the applied surface 
current. 

For a large number of cases, we have performed both 
calculations and found results in very good agreement. 
However, the first approach, while computationally much 
less expensive, shows, due to the open boundary condi- 
tions used in the DMRG runs, much more pronounced 
finite-size effects: SF/DDW correlations are suppressed, 
not enhanced by the presence of open boundary condi- 
tions. Long-ranged or quasi-long-ranged order is thus 
often much harder to detect. Figs. [S] and show 
the current-current correlations for a point in parame- 
ter space which exhibits SF/DDW order both for 8 and 
12 percent doping and various ladder lengths. In the first 
case, correlations are weakly suppressed by the edges and 
order is demonstrated in agreement with the second ap- 
proach; in the latter, correlations arc strongly suppressed 
by the edges and grow strongly in size with ladder length; 
even for 200 rungs, it would be impossible to detect long- 
ranged order which is clearly visible from 50 rungs up- 
wards in the second approach. For parameter sets where 
accessible system sizes allow to see long-ranged correla- 
tions also in the rung-rung current correlations, the val- 
ues of the rung currents that can be extracted are in 
very good agreement with the values found by imposing 
an infinitesimal edge current. 

In the second, computationally much harder, ap- 
proach, we also observe that on very short ladders right at 
the CDW-SF/DDW phase boundary the open boundary 
conditions strongly favor a checkerboard charge density 
wave over the SF/DDW currents (Fig. UJ, but the mag- 
nitude of the plaqucttc currents stabilizes very quickly 
with increasing system length, such that we have pre- 
ferred this approach in our numerical calculations. In 
Refs. I4fj|l4ll exact duality relations were used to ascer- 
tain that the phase boundary between the SF/DDW and 
CDW phases lies at U - V± + (3/4) J± = 0. It is re- 
markable that, as evident in Fig. [7] strong currents are 
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FIG. 5: Rung current-current correlations C(r) calculated for 
ladders of length 50 (solid circles), 100 (open circles), 150 
(diamonds) in the SF/DDW phase for U = -0.4, t = 1, 
V± = 0.9, J± = 2, and doping 5 = 0.08. Jy = Vj| = in this 
and all other figures in this section. 
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FIG. 6: Rung current-current correlations C(r) as in Figure 
but doping 5 = 0.12. 



found to persist right at the phase boundary separating 
the SF/DDW and CDW phases. In this particular calcu- 
lation, we have taken the edge current as small as 1CP 8 < 
without observing any change in the results. 

The SF/DDW phase can be clearly distinguished from 
the neighboring CDW and doped D-Mott phases^ In 
both the CDW (Fig. 0) and doped D-Mott (Fig. 
phase, the edge current induces rung currents j± (x) that 
die out exponentially fast, while in the SF/DDW phase it 
stabilizes a rung current orders of magnitude larger (Fig. 
ITU1) . In the SF/DDW phase the staggered rung current 



FIG. 7: Rung current at the phase boundary between 
SF/DDW and CDW order. The edge current is O.OOOlt and 
U = 0.25, t = V± = 1, Jx = 1, and doping 8 = 0.04 for ladder 
lengths 50 (open circles) and 100 (solid circles). 



j s (x) (defined in Eq. (|3.16|) ) oscillates around zero with 
a characteristic wavelength 2/6 (this wavelength tends to 
decrease slightly with decreasing distance to the bound- 
ary). Thus the rung currents have anti-phase domain 
walls separated by a characteristic distance 1/5. At these 
domain walls (or "discommensurations" ) the rung cur- 
rent has the same sign on adjacent rungs (H — h or ) 

instead of alternating between rungs as in the ( "commen- 
surate") regions between domain walls. Furthermore, us- 
ing that 



2kp — kpi + kp2 = nn = 7r(l — 5), 



(4.5) 



which is the generalization of Eq. 1)3.1(1 away from half- 
filling, one observes that the rung current oscillations 
may be associated with the fundamental wavevector 2k p 
as these oscillations vary approximately as cos2kpx = 
(— l) x cosirSx (this equality holds since the rung index x 
is an integer). Moreover, the currents on the top and bot- 
tom leg are always opposite in sign. Consequently, the 
ordering wavevector of the currents is (2kp, n), which re- 
duces to (-7T, 7r) at half-filling, as expected. 

A modulation of the rung charge density coexists with 
the orbital currents in the SF/DDW phase (Fig. [TUj l. 
The characteristic wavelength of this modulation is 1/6 
(again, the wavelength decreases slightly as the distance 
to the boundary decreases), which is just the average 
distance L/(A/ r ? l /2) between pairs of doped holes. Thus 
there are two doped holes per wavelength, one per leg; in 
a cartoon picture the doped hole pairs simply distribute 
themselves equidistantly along the ladder, each hole pair 
occupying one rung. These hole pairs tend to stay in the 
vicinity of the rung current minima (the domain walls of 
the current pattern). The ordering wavevector for this 
charge density modulation is (4/cf,0) because its charac- 
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FIG. 8: Rung current (solid circles, left j/-axis, units of t) 
and electronic densities on top and bottom ladder leg (up 
and down triangles, right j/-axis) for 100 rung-ladder in the 
CDW phase and edge current 0.000 H (17 = 0.25, t = Vi. = l, 
J±_ = 0.8, and doping 5 = 0.04). 



teristic periodicity in the direction parallel to the legs is 
the same as that of cos4fcp2; = cos27t&e, and the density 
is the same on the top and bottom site of a rung. Clearly, 
this charge density modulation disappears at half-filling. 

As seen in Fig. E3 the (4fcp , 0) modulation of the charge 
density exists also in the doped D-Mott phase. The na- 
ture of the charge density pattern in the CDW phase 
(Fig. |5J is less obvious at first sight; however, as will be- 
come apparent in Sec. it is the sum of a large compo- 
nent with wavevector (2fcp,7r) and a smaller component 
with wavevector (4/cf,0). Thus the (4fcp,0) charge den- 
sity modulation is present in all three phases appearing 
in the parameter region considered here. Such (4/cp,0) 
charge density modulations in two-leg ladders have re- 
cently been discussed by White et al.^ we will analyze 
these modulations in more detail in Sec. ]V\ 

We have also calculated the DSC correlations in the 
SF/DDW and doped D-Mott phases (Fig. [TTJ>. In the 
SF/DDW phase the DSC correlations decay exponen- 
tially. As expected, the DSC correlations fall off more 
slowly in the doped D-Mott phase. Close to the phase 
boundary to the SF/DDW phase the correlations appear 
still to be exponential, however, while further into the 
doped D-Mott phase the edge effects make it difficult to 
say whether for longer system sizes the correlations might 
turn out to be algebraic. 

Most of our DMRG calculations are for particle densi- 
ties close to half-filling. As DMRG is a canonical ensem- 
ble method, particle numbers are fixed at each growth 
step and there is some arbitrariness in inserting the par- 
ticles and holes to maintain essentially constant particle 
density during system growth. Various insertion schemes 
have yielded essentially identical results. To obtain such 



FIG. 9: Rung current (solid circles, left j/-axis, units of t) 
and electronic densities on top and bottom ladder leg (up 
and down triangles, right y-axis) for 100 rung-ladder in the 
doped D-Mott phase and edge current O.OOOli (U = 0.25, 
t = V± — 1, J± = 1.7, and doping S = 0.04). Top and bottom 
densities coincide. Note that for these model parameters there 
is SF/DDW order at half-filling (see Fig. |Ua)); the phase 
boundary between the D-Mott and SF/DDW phase moves 
downwards when the system is doped. 



well-converged results, it is crucial to apply several runs 
(sweeps) of the finite-system DMRG algorithm until con- 
vergence is observed numerically. We have typically per- 
formed of the order of 10 to 20 sweeps and observed con- 
vergence after typically 5 to 7 sweeps. 

Up to 800 states were kept for DMRG runs; follow- 
ing standard practice, the number of states kept was ini- 
tially chosen to be much smaller and increased during 
sweeps. For systems up to about 100 rungs, final results 
were independent of the way the number of states was 
augmented. For systems well beyond 100 rungs we have 
observed that in the SF/DDW phase, DMRG runs that 
started with a very low number of states kept, converged 
to a result where the current amplitudes were somewhat 
suppressed with respect to runs that started with high 
precision and yielded results in perfect agreement with 
those found for shorter systems. The low precision cal- 
culation introduces a phase slip it at the center (Fig. I12f> . 

Do we really have true long-range order in the orbital 
currents? If there is not true long-range order, one might 
speculate that the introduction of a phase slip (of the 
type seen in Fig. I12f) at low energetic cost is the way how 
long-ranged correlations are lost in this system. How- 
ever, we did not observe, for systems up to 200 rungs, 
any significant decrease of the magnitude of the trig- 
gered current pattern. Considering Fig. 1131 one actually 
sees that maximum current amplitudes grow from system 
size 50 to 100, in line with the SF/DDW suppression by 
edges reported above, to stay constant after that. One 
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FIG. 10: Rung current (solid circles, left y-axis, units of t) 
and electronic densities on top and bottom ladder leg (up 
and down triangles, right j/-axis) for 100 rung-ladder in the 
SF/DDW phase and edge current O.OOOlt (U = 0.25, t = 
V± = 1, Jx = 1.5, and doping S = 0.04). Top and bottom 
densities coincide. 
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FIG. 12: Staggered rung current for 200 rung-ladder in the 
SF/DDW phase and edge current O.OOOlt (U = 0.25, t = 
V± = 1, J± = 1.5, and doping S = 0.04). Final number of 
states 77i = 400 was reached starting from m = 100 (open 
diamonds) and m = 400 (solid diamonds) . 
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FIG. 11: Log-log plot of DSC correlations for a 100-rung 
ladder with U = 0.25, t = V± = 1 and doping S = 0.04, 
for different values of J± (from bottom to top): 1.3 (solid 
squares), 1.5 (open diamonds) , 1.7 (solid diamonds), 2.0 (open 
squares), and 2.6 (solid circles). The first two values are in the 
SF/DDW phase, the last three in the doped D-Mott phase. 
The inset shows a lin-log plot of the same data. 



might envisage that the open boundary conditions lock in 
the currents, leaving us with seeming long-ranged order. 
The repeated observation that open boundary conditions 
seem to disfavor SF/DDW does not support this point 
of view. It might also be speculated that our filling fac- 
tors, 4, 8 and 12 percent away from half-filling, allow the 
formation of commensurate current patterns which may 
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FIG. 13: Rung current in the SF/DDW phase and edge cur- 
rent O.OOOlt (U = 0.25, t = Vx = 1, J± = 1.5, and doping 
S = 0.04) for ladder lengths 100 (solid squares), 150 (circles), 
200 (diamonds). 



show long-ranged order as opposed to the generic case 
of incommensurate filling. Introducing 11 holes on a 192 
rung ladder to model an approximation to a generic in- 
commensurate filling on a finite system, however, we do 
not observe current decay either (Fig. I14|) . 

The SF/DDW phase survives up to quite substantial 
dopings away from half-filling. Fig. H5l shows current and 
density oscillations in the SF/DDW phase of a 200-rung 
ladder for doping S = 0.08. For the parameter sets stud- 
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FIG. 14: Rung current (open squares, left y-axis, units of t) 
and electronic densities (solid circles, right y-axis) for 192- 
rung ladder in the SF/DDW phase (U = -0.4, t = 1, V± = 
0.9, J± = 2), edge current 0.000 U and 11 holes (i.e., 5 ss 
0.0286, 1/8 « 34.9). The disturbed density pattern and slight 
current reduction on the far end of the ladder is due to the 
presence to one unpaired hole. 
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FIG. 15: Rung current (lines, left j/-axis, units of t) and elec- 
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ladder in the SF/DDW phase for 5 — 0.08 and edge current 
O.OOOlt (U = 0.25, t = V± = 1, Jx = 1.5). 



ied, the phase transition where SF/DDW order is finally 
suppressed to zero seems to occur roughly between 10 
and 20 percent doping; a more precise estimate would re- 
quire a more systematic study of the doping dependence 
of the SF/DDW phase as a function of the interaction 
parameters. 



V. EFFECTIVE FIELD THEORY 
DESCRIPTION OF THE DOPED LADDER 

In this section we try to understand the DMRG results 
in Sec. IIVI using weak-coupling RG/bosonization argu- 
ments. We will be particularly interested in what these 
arguments predict regarding the possibility of having true 
long-range order for the considered doping levels. 



A. Physical observables 

We first present the set of relevant physical observ- 
ables and their bosonized expressions for the case of an 
arbitrary filling for which the Fermi energy crosses both 
bands. This includes the half-filled case discussed in Sec. 
IIIII as well as the doping levels studied numerically in 
Sec. HV| 

The rung current operator is 

3±(x) = 32k P [ cos(2fci?a; + <j> +p ) cos 4> +(T cos 6- p cos 
+ sin(2fci?a; + 4>+p) sin 4>+ a sin £9_ p sin 0- a ~\ . (5-1) 

Furthermore, the density operator ng(x) on leg I is 

IT OX 

± n2k F [ cos(2kpx + 4>+ p ) cos <j) +a sin 9- p cos 6_ a 

— sin(2kpx + <fi+ p ) sin <fi +a cos #_ p sin (9_ CT ] 

+ n^kp cos(4kpx + 2(f> +p ) cos2</> +(T . (5.2) 

The quantity 2kp was defined in Eq. (|4.5|) . The nonuni- 
versal coefficients j2k F , ^2fe F and depend on the 

short-distance cutoff of the theory. 

The Akp term in ng{x) should be particularly noted. 
It does not come out of a naive calculation of the density 
operator using the bosonization formula 1)3. 7J) : however, 
a more general, phcnomcnological reasoning shows that 
such higher harmonics terms are generally expected*^^ 
The form of the 4&f term in nt{x) was deduced by White 
et aZ4£ They considered the AkF term in the correlation 
function for the square of the density operator, as result- 
ing from the product of the 2kp terms in the density 
operator calculated from the bosonization formula, and 
then (implicitly) argued that a similar Akp contribution 
to that correlation function would be expected to come 
from the product of the constant term and a AkF term 
in the density operator. Higher harmonics than 2kF in 
jj_(x) and AkF in ni(x) have been neglected in the ex- 
pressions above as they are expected to give at most only 
minor quantitative corrections to the terms already in- 
cluded. 

The density operator above is more complicated than 
the one discussed for the half-filled CDW phase in Sec. 
IIII Dl It is therefore instructive to see how the expecta- 
tion value of ni{x) reduces to the correct form at half- 
filling. In that case, both <p +p and <j>+„ are locked in 
all four phases (see Table QJ. This gives (d x <p+ P ) = 



12 



and (cos20 +(T ) ^ 0; furthermore, from Akp = 2ir it fol- 
lows that {cos{^kpx -\- 2(f) +p )) = (cos2(/> +p ) ^ 0. Thus 
the expectation value of the 4fcp term is independent of 
x at half-filling. Consequently, the only way to obtain 
the correct average density of one electron per site is to 
have riik F = in this case. Hence the expectation value 
of the deviation of the density from its average value, 
Sni(x) = ni(x) — 1, reduces to the expectation value of 
the 2kp term, which is only nonzero in the CDW phase, 
in agreement with the simplified discussion of the density 
operator in Sec. IIII Dl 

Finally, we note that since the DSC and SSC operators 
have zero momentum, their bosonized expressions are the 
same as at half- filling; see Eqs. (|3.18jl - H3.20jl . (We will 
not consider SSC order in the following, since the doped 
S-Mott phase was not studied in Sec. II VI ) 

B. True long-range order scenario 

It was argued in Sec. IIVI that the DMRG results 
are consistent with true long-range orbital current and 
charge density wave order. From the bosonization point 
of view, to have true long-range order it is necessary 
that certain bosonic fields become locked to appropriate 
values. In this subsection we show that if this locking 
occurs, the resulting expectation values and correlation 
functions of the physical observables of interest, as cal- 
culated from the bosonized expressions in Sec. IV Al arc 
in good qualitative agreement with the DMRG results. 
The assumption that the locking occurs is easily justified 
for all bosonic fields in question except the symmetric 
charge field <fi+ p - The conditions for locking this field 
will be investigated in detail in Sec. IV CI 

Let us first consider the SF/DDW phase. True long- 
range order in the currents, (j±(x)} ^ 0, requires that 
the four fields <fi+ p , (f>+ a , and 9^ a become locked 

to appropriate values as determined by Eq. 15.111 . 
The momentum-conserving interaction (|3.13|) , which con- 
tains cosines of the three latter fields, is present also 
away from half-filling, and weak-coupling one-loop RG 
calculations^iiSSiSLS show that it can cause 4>+ a , Q-p, 
and 9 - a (or (t>- a ) to lock also in the doped case, at least 
if the doping is not too large. This picture is expected to 
hold also for stronger interactions. Let us now assume 
that also the symmetric charge field </>+ p can become 
locked for the dopings considered in Sec. IIVI There 
is then a finite energy gap to excitations in </>+ p . Using 
the field-locking pattern for </>+<t, 9- p , and in Ta- 
ble[J and evaluating the expectation value of the current 
operator scmiclassically^ we find 

(j±(x)) ~ jo cos(2fc f .T + (</> +p )), (5.3) 

with jo = j2k F (cos 4>+o- cos 9-p cos 9- a ) ^ 0. Since +p 
and </>+ CT are locked, we also have (with nt{x) = rib(x) = 
n(x) in the SF/DDW phase) 

(n(x)) = (1 - 5) + n cos(4k F x + 2{<j> +p }), (5.4) 



with rig = n4k F (cos 2(j) +a ) ^ 0. This shows that (2fcj?,7r)- 
SF/DDW order and (4fc F ,0)-CDW order coexist in the 
SF/DDW phase. Furthermore, if the non-universal coef- 
ficient riikp is positive, the minima of the electron density 
always occur at the zeros (anti-phase domain walls) of the 
current. These results are in agreement with the findings 
in Sec. IIVI The different solutions for (<fr+ p ) give rise to 
degenerate ground states which are related to each other 
by translations by a lattice vector (in the DMRG calcu- 
lations, this degeneracy is of course lifted by the open 
boundary conditions, leaving only a two-fold degeneracy 
due to the breaking of time reversal symmetry). Fig. 1161 
shows a fit of these analytical results for the SF/DDW 
phase to the DMRG results for two different dopings. 
The agreement is quite good. If the true long-range or- 
der scenario is in fact realized, we expect that the minor 
differences between the numerical and analytical results 
could be removed, at least in principle, by improving the 
latter by taking into account boundary and finite-size ef- 
fects, corrections due to the presence of higher harmonics, 
and by going beyond the semiclassical approximations 
used in evaluating the expectation values. 

By using similar arguments as for the SF/DDW phase, 
one also predicts, again in agreement with DMRG results, 
that there is coexistence of (2kp, 7r)-CDW and (4fcp,0)- 
CDW order in the CDW phase, while in the doped D- 
Mott phase there is "only" (4fc F ,0)-CDW order. The 
reason for the ubiquity of (4fcp,0)-CDW order is that 
this order parameter only contains the fields <p+ p and 
4>+a , both of which are locked in all phases in this true 
long-range order scenario. 

This scenario furthermore predicts that the asymptotic 
decay of DSC correlations is exponential in all phases. 
The DSC correlation length is expected to be smallest 
in the CDW phase because then the field-locking differs 
from what the DSC operator wants in both the —p, —a, 
and +p sectors. In the SF/DDW phase the correlation 
length should be larger because there is no longer dis- 
agreement in the — p sector. In the doped D-Mott phase 
the correlation length is larger still because the —a sector 
also matches up, the only exponential decay left coming 
from the +p sector. These predictions for the qualitative 
behavior of the DSC correlations are in good agreement 
with the DMRG results (with the qualifier that in the 
case of the doped D-Mott phase, the nature of the asymp- 
totic decay of DSC correlations is difficult to deduce from 
the DMRG data; see the discussion in Sec. IIV() . 



C. True long-range order and Umklapp 
interactions 

The crucial assumption which enabled the true long- 
range order scenario just discussed to emerge from the 
bosonization formulas was that the symmetric charge 
field (p+p can become locked for the dopings and inter- 
action parameters considered in Sec. IIVI In this subsec- 
tion we use symmetry arguments to derive the form of 
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Rung 



FIG. 16: (Color online) (a) Staggered rung current j a (x) (red 
curves, left y-axis) and electron density per site n(x) (blue 
curves, right y-axis) as a function of rung index x in the 
SF/DDW phase of a 200-rung ladder with 8 = 0.02. Shown 
are DMRG results (diamonds and circles for currents and den- 
sities, respectively, connected by dashed lines to guide the eye) 
and analytical fits from the true long-range order scenario (full 
lines). The model parameters for the DMRG calculation are 
t = V± = 1, U = 0.25, J± = 1.5 and Jy = Vj| = 0. The values 
of the fitting parameters are |jo| ~ 0.61 and no ~ 0.0135. (b) 
Plots for the same model parameters as above, except that 
the doping S — 0.04 and the values of the fitting parameters 
are \j \ ~ 0.52 and n « 0.0132. 



the terms that, under the right circumstances, are able 
to lock </>+p. It will be shown that these terms are Umk- 
lapp interactions, and that they are allowed by symme- 
try in the low-energy effective Hamiltonian only when 
the average number of electrons per site, n, is a ratio- 
nal number p/q. The further analysis suggests that the 
Umklapp interactions are only able to lock </> +p for rel- 
atively small values of q. Therefore this weak-coupling 
RG/bosonization analysis appears not to support a true 
long-range order scenario for the dopings considered in 



Sec. IIVI which correspond to quite large (> 25) values 
of q. 

If the low-energy effective Hamiltonian, denoted by 
.ffeff, is invariant under global translations in 0+ p , 

<pj rP {x) — y <pj rP {x) + a (a arbitrary real), (5-5) 

locking </>_|_p implies that this continuous symmetry must 
be spontaneously broken, which is impossible in one di- 
mension. Therefore, to lock <p +p it is necessary that 
there be terms in H c g that break this continuous sym- 
metry. Furthermore, the Hamiltonian can only contain 
gauge-invariant terms. The bosonic fields themselves are 
not gauge-invariant, only derivatives and exponentials of 
them are. Of these, only exponentials of <fi+ p are not 
invariant under i|5.5[) . as they break the continuous sym- 
metry down to a discrete one. Thus one is led to consider 
terms of the form 

H Q =F Q + F q , F q ^ Jdx Aq (x) e^+oW (5.6) 

where Q ^ and Aq (x) is an operator that we define to 
be invariant under 1)5.5(1 . We may take Q > without 
loss of generality. 

We will use the fact that Hq must be invariant under 
certain symmetry operations (reflections and translations 
in real space) to deduce under what conditions Hq can 
appear in H c ff, and what the allowed values of Q are in 
that case. However, to do this it turns out that we will 
also need to know how Hq transforms under translations 
in 4>- p . Since the conjugate field 8^ p is locked in all 
phases, 0_ p is disordered, so that an operator containing 
terms with exponentials of 0_ p will be irrelevant in the 
RG sense. Thus, if Hq is to have a chance of causing true 
long-range order, we must assume that Aq{x) does not 
contain such terms, which implies that Hq is invariant 
under global continuous translations in (f>- p , 

4>-p — > 4>- p + (3 (/? arbitrary real). (5.7) 

Let us first consider reflection symmetry. Hq must be 
invariant under the transformation 

Ct,x,s * Cb,x,s? Cb.x.s * Ct.x.s (^-^) 

which interchanges the "top" and "bottom" leg. By 
regarding the two sites on a rung as forming a one- 
dimensional lattice in the transverse (y) direction with 
periodic boundary conditions, we see that (|5.8|l is a trans- 
lation by one lattice spacing in this two-site lattice, and 
is therefore effected by the unitary operator e~ lVy where 
V y is the momentum operator in the y direction. The 
allowed transverse momenta arc k y = and 7T, associ- 
ated with the bonding and antibonding bands 2 and 1, 
respectively. Thus 

' p y= 7: J2 c ix S c ixs=TrJ2 Afpis - ( 5J ) 

x , s P, s 
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One finds that e %Vy transforms the continuum field op- 
erators as 

1ppis{x) — > —^pu{x), lpP2s(x) — > 1pP2s(x). (5.10) 

Thus the symmetry under leg interchange implies that 
any term in the Hamiltonian must contain an even num- 
ber of both 1- and 2-operators of this type (the latter 
due to conservation of the total number of particles). 
Bosonization gives V y = J dx (d x <fi+ P — d x <f>- p ) which 
shows that e~ lVy generates translations in 0± p . We would 
however like to find an operator that transforms the field 
operators like l|5.10[) while at the same time generating 
translations in <f>± p , as that would give information about 
Q. Such an operator is e~ llz , where R is obtained by 
modifying V y in (|5.9() by a factor of P, which changes 
<p — > 9 in the bosonic representation: 

n = nJ2 p Mpu = J dx (d x e +p -d x 0- p ). (5.11) 

P,s 

As desired, e~ lli leaves all bosonic fields invariant except 
cf>±p which transform as 



7T, <f>. 



+p 



(5.12) 



Thus invariancc of Hq under this transformation implies, 
when we use that Hq must be invariant under (|5.7|) sep- 
arately, that Q must be an even integer. 

Next, consider translation symmetry^ Hq must be 
invariant under the transformation 



ce,x+i,s 



(5.13) 



which translates the system by one lattice spacing in the 
x direction. This transformation is generated by the op- 
erator e -1 ^ where V x is the momentum operator in the 
x direction, given by 



V. r 



— I dx 

7T 



E 

PAs 



= 2J k n k>. 

kXs 

[(fc.Fi + k F2 )d x 9 +p 



Pk F xNpx s 

{k F2 - k F1 )d x 0. 



(5.14) 



Clearly e i ' Px leaves all bosonic fields unchanged except 
(j)± p which transform as 



'-p 



'-p 



k F i 
k F \ 



k F2 , 
k F2 



'+p 



(5.15) 
(5.16) 



where we used (|4.5ll to express <|5.16ll in terms of the 
average number of electrons per site, n. Again invoking 
the invariance of Hq under (|5.7|) . it follows that (|5.1t)|) 
requires Qn to be an even integer. This is impossible 
if n is an irrational number, implying that Hq is not 
allowed in H c ^ in this case. Therefore true long-range 
order is excluded for incommensurate fillings. If instead 
n is rational, i.e., 



n = p/q 



(5.17) 



where p and q are coprimes (i.e., integers with no common 
divisors), the values of Q that make both Q and Qn even 
integers are 



Q = 2rq, r = 1,2, 



(5.18) 



Therefore the terms H 2rq are allowed in H c r in this case. 

We now consider the physical interpretation of these 
terms. As they are not invariant under (|5.5|l . they do not 
commute with the generator of these translations, given 
by G 4> +P = I dx n 0+ P = (A/i - AAr)/2, where IT 0+p = 
—dxO+p/n is the conjugate momentum of </>+ p , and 

N P = J2Afpx s = - [dx d x {^ +p + P6 +p ) (5.19) 

As 71 

is the number operator for right-/left-moving electrons. 
Because A2 rq (x) is invariant under independent global 
continuous translations of <p +p and 6 +p (the latter due 
to conservation of total particle number), we have 
[Af P , A 2rq (x)} = 0. This gives 



[Afp, F< 



2rq\ 



-2rqP F 2 



rq • 



(5.20) 



which shows that F 2rq and F\ rq are raising/lowering op- 
erators (by 2rq quanta) for Afp: F 2rq scatters 2rq elec- 
trons from right to left while F\ rq scatters them the other 
way. Furthermore, the commutator 



[V X ,F 



2rq\ 



-4k F rq F 2 . 



rq 



(5.21) 



shows that the associated momentum transfer is ^fAkprq 
for F 2rq and F^, respectively, which reduces to ^f2nrp 
for n = p/q. Thus H 2rq is an Umklapp interaction. Note 
that the invariance of Hq under l|5.7|l was used in the 
derivation of (|5.21|) . 

Of all the Umklapp interactions H 2rq , r = 1,2, . . ., that 
arc allowed in H c h for commensurate density n = p/q, 
the operator H 2q corresponding to r = 1 is most relevant 
in the RG sense, and is therefore the one we will focus 
on in the rest of the discussion. H 2q scatters 2q electrons 
from the left to the right side of the Fermi surface (and 
vice versa), with associated momentum transfer ±27rp£ 
The interaction H 2q with q = 1 (which is allowed in H c g 
at half-filling) will be referred to as the basic-Umklapp 
interaction. The interactions H 2q with q > 1 will be 
referred to as multiple- Umklapp interactions. 

The derivation presented here shows that when the 
density n is moved away from the commensurate value 
p/q, the term H 2q is no longer allowed in H a f[. Another 
way of understanding this is from the fact that when n 
moves away from p/q, A 2q (x) acquires factors that oscil- 
late in space (this is intimately related to the fact that the 
momentum Akpq transferred by H 2q now deviates from 
the reciprocal lattice vector 2-Kp), and these make H 2q 
average to zero on sufficiently long length scales, thereby 
rendering it an irrelevant operator. For example, for dop- 
ing S away from half-filling the basic-Umklapp interaction 
is proportional to cos(20 +(9 + 4fc^a;) = cos(2<j> +p + 2ir5x). 
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Thus the cosine oscillates with x and averages to zero 
over length scales much bigger than 1/5. 

Since the —p and ±cr sectors are gapped out, the 
parts of the Hamiltonian containing fields belonging to 
these sectors can be replaced by their expectation val- 
ues. The low-energy effective Hamiltonian is then de- 
scribed in terms of the fields in the remaining +p sector. 
H e g has a quadratic part given by the Luttinger liquid 
HamiltonianS 



H Lh = I dx 



2tt 



K+p(dx9+p) 2 + -^—(d x 4>+p) 



, (5.22) 



where v +p and K +p are the effective velocity and Lut- 
tinger parameters, respectively. i?LL describes a crit- 
ical (Luttinger liquid) fixed point with gapless excita- 
tions in (f>+ p . For commensurate density n = p/q, H e g 
will also contain H 2q if the latter is a relevant operator 
with respect to Hll- In that case H 2q will lock (j) +p 
in one of the minima of its potential, thereby opening 
up a gap in (j>+ p , causing true long-range order. Assum- 
ing that (A 2q (x)) 7^ for the field-locking pattern in 
the gapped sectors (if (A 2q (x)) should happen to van- 
ish, true long-range order is ruled out), the only decay of 
the two-point function of H 2q = J dx 7i 2q (x) comes from 
the exp (±2i<70 +p ) factors, giving the asymptotic decay 
{H 2q {x)H 2q {0)) ~ \x\~ 2d ^, where 



d 2q = q 2 K. 



-p 



(5.23) 



is the scaling dimension of H 2q . If d 2q < 2, H 2q is pertur- 
batively relevant. Thus the critical value of K +p below 
which H 2q is expected to open up a gap and cause true 
long-range order is from this simple argument estimated 
to be 



Kp = Vi 2 



(5.24) 



The dopings considered in Sec. IIVI correspond to very 
large values of q, and consequently very small values of 
K+ p (see Table [Hj. For weak repulsive interactions K +p 
is expected to be slightly less than 1 (K +p = 1 in the 
noninteracting model). Therefore H 2q is strongly irrele- 
vant for weak interactions in the doped ladder, and conse- 
quently true long-range order is not possible in that case. 
For larger interactions a definite statement can no longer 
be made, as K +p is not known as a function of U, Vj_, J± 
etc. On the other hand, true long-range order would cer- 
tainly be very surprising. Although one generally expects 
K+p to decrease with increasing interaction strength, the 
values of K +p for the nonzero dopings in Table [H] are, to 
the best of our knowledge, orders of magnitude smaller 
than the smallest known values of K +p in the t-J and 
Hubbard ladders^ which furthermore are obtained only 
for very strong interactions, much stronger than the ones 
used in Sec. IIVI Correspondingly, only for low commen- 
surabilities (i.e., small values of q) has it been found that 
K +p can become as small as K +p for sufficiently strong 
interactions in theoretical ladder models, with true long- 
range order as a result. 



S(%) 


P 


q 


K +P 





1 


l 


2 


2 


49 


50 


0.0008 


2.86 


373 


384 


0.0000136 


4, 8, 12 


24, 23, 22 


25 


0.0032 



TABLE II: The parameters p, q, and K\ p — 2/q 2 evaluated 
for the rational dopings 5 = 1 — p/q considered in Sec. IIVI 
The values of these parameters at half-filling are shown for 
comparison. 



We end this subsection with some remarks on how the 
discussion presented here relates to previous work. It 
is already known in the literature that RG/bosonization 
arguments predict that having true long-range order in 
the two-leg ladder requires relevant Umklapp interac- 
tions and that Umklapp interactions are usually irrel- 
evant away from half-filling. These conclusions, and 
arguments leading to them, have been invoked either 
explicitly or implicitly in many previous studies of in- 
teracting electrons on a two- leg ladder (see, e.g., Refs. 

The discus- 
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sion given here based on symmetry considerations, which 
is a generalization to the two-leg ladder of previous dis- 
cussions for a single chain allowed for a more sys- 
tematic and general derivation of these results. In par- 
ticular, the derivation of Eq. (|5.23|) enabled us to give a 
semi-quantitative discussion of the question of how irrel- 
evant (or relevant) the most relevant allowed Umklapp 
operator is for an arbitrary rational filling. 



D. Quasi long-range order scenario 

The weak-coupling analysis in the previous subsection 
suggests that Umklapp interactions are unlikely to be 
relevant for the dopings considered in Sec. IIVI and that 
as a consequence the effective theory at sufficiently low 
energies and long wavelengths is given by the Luttinger 
liquid Hamiltonian in Eq. i|5.22[l . with gapless excita- 
tions in 4>+ p . This will be referred to as the quasi long- 
range order scenario, as it implies that for sufficiently 
large distances correlation functions of exponentials of 
<j)+p decay as power laws with fT +p -dependent exponents 
(the same features hold for 9 +p ). In this subsection we 
summarize some consequences of this scenario, both for 
infinite and finite ladders. From a numerical compari- 
son with the Luttinger-liquid predictions for finite lad- 
ders, we find that the DMRG results are not consistent 
with the quasi long-range order scenario. This analysis 
complements the discussions in Sees. IIVI and IV Bl which 
concluded that the DMRG results are consistent with a 
true long-range order scenario. 

We start by summarizing the Luttinger-liquid predic- 
tions for the correlation functions for the charge density, 
orbital current and DSC pairing operators in the doped 
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Correlation 
function 


Phase 


SF/DDW 


CDW 


Doped D-Mott 


(2fc F ,7r)-SF/DDW 
(2fc F ,7r)-CDW 

DSC 
(4fc F ,0)-CDW 


K+p/ygL (const) 
exp (exp) 
exp (exp) 
2K+p (const) 


exp (exp) 
A' + p/2fiZ. (const) 
exp (exp) 
2K+P (const) 


exp (exp) 
exp (exp) 
1/(2A +P ) 31 . 37 . 45 i 80 (exp) 
2K +P 3±^2. (const) 



TABLE III: Bosonization predictions for the asymptotic decay of two-point correlation functions in the two-leg ladder. In each 
cell, the first entry is the decay when 4>+p is gapless (quasi long-range order scenario); the second entry (in parenthesis) is the 
decay when cj>+p is gapped (true long-range order scenario). Power-law decay is indicated by the decay exponent, 'const' means 
no decay, and 'exp' means exponential decay. 



SF/DDW, CDW and D-Mott phases in an infinitely long 
ladder. The nature of the decay of these correlation func- 
tions is easily obtained from the bosonized expressions in 
Sec. IV Al the field-locking patterns for the ±er and — p 
sectors (Table HJ, an d the correlation functions of expo- 
nentials of 4> +p and Q+p. The results are summarized in 
Table ITTT1 A few remarks are in order. The 4&f density 
correlations show the same qualitative behavior across 
the entire phase diagram. In the SF/DDW and CDW 
phases, the SF/DDW and 2/cf-CDW correlations, re- 
spectively, dominate over the 4fcp-CDW correlations. In 
the doped D-Mott phase, the DSC correlations are dom- 
inant for K +p > 1/2 while the 4&F-CDW correlations 
dominate for K +p < l/2. 31 i 45 i 80 Schukja has argued that 
K +p — > 1 as 5 — > + in the doped D-Mott phase. There- 
fore one expects DSC correlations to be dominant in the 
doped D-Mott phase sufficiently close to half-filling. (For 
easy comparison the nature of the asymptotic decay when 
4>+ p is gapped has been indicated in parenthesis in Ta- 
bic ^Oj The cases we have commented on here are those 
for which the true and quasi long-range order scenarios 
predict qualitatively different decays.) 

The expressions for the power-law exponents in Table 
IIIII spur the question of whether the doped ladder might 
be described by a quasi long-range order scenario with a 
very small K +p . Unfortunately, the strong finite-size ef- 
fects in the DMRG results severely complicate a compar- 
ison with the Luttingcr-liquid predictions for an infinite 
ladder. It is much preferred to compare the DMRG re- 
sults with analytical predictions for a ladder of finite size 
and with open boundary conditions. White et alm^ have 
recently discussed such a ladder that is in a Luttinger- 
liquid phase of the type described above, i.e., in which 
the low-energy effective theory is given by the gapless 
Luttinger- liquid Hamiltonian i?LL in the +p sector and 
the three other sectors are gapped. It was shown that 
in such a system the boundaries will act as impurities 
and induce generalized Friedel oscillations in the charge 
density which decay away from the boundaries. These 
Friedel oscillations have fundamental wavevector 4fcp and 
thus period 1/5. White et al. furthermore showed that 
the density amplitude in the middle of the ladder should 
scale with the ladder length L as L~ +>•. To investigate 
whether the 4fcp charge density oscillations found in Sec. 



HVI might be interpreted as generalized Friedel oscillations 
in such a Luttinger liquid, one can plot the mid-ladder 
density amplitude as a function of the ladder length in a 
log-log plot and try to fit the results to a straight line; 
from its slope one can extract a value for K +p m^ 

We have attempted this kind of fit of the DMRG results 
for two different cases, both at 4 percent doping: In the 
SF/DDW phase {t = V± = 1, U = 0.25 and J ± = 1.5), 
and in the doped D-Mott phase (same parameters ex- 
cept J± = 1.8). In both cases the mid-ladder density 
amplitude is found not to decay with ladder length for 
sufficiently large ladders, so that the fit gives K +p = 0. 
This result does not support an interpretation of the 
density oscillations as generalized Friedel oscillations in 
a Luttinger liquid; instead it is consistent with a true 
long-range order scenario. Even if one were to take into 
account corrections to Hll due to irrelevant operators 
(such as the basic-Umklapp interaction discussed more 
below), one would still expect that the mid-ladder density 
amplitude should decrease monotonically towards zero 
with increasing ladder length, instead of approaching a 
nonzero constant. Thus if the DMRG results reflect the 
true behavior of the system the quasi long-range order 
scenario appears to be ruled out. 



E. Neglecting quantum fluctuations: A classical 
description of the doped ladder 

The fact that the basic-Umklapp interaction is rele- 
vant in the half-filled ladder enabled true long-range or- 
der in that case. The analysis in Sec. IV CI suggests that 
in the doped ladder the basic-Umklapp interaction is ir- 
relevant: although it can cause significant ordering ten- 
dencies on length scales < 0(1/6), it can be neglected 
on length scales 3> 1/5 The Luttingcr-liquid behav- 
ior discussed in Sec. IV Dl with quasi long-range ordered 
correlation functions, would thus be expected to mani- 
fest itself in this asymptotic regime^ One would think 
that this regime should be accessible numerically since 
the dopings considered in Sec. IIVI correspond to values 
of 1 / 5 that are considerably smaller than the largest sys- 
tem sizes studied there. In spite of these expectations, 
we have seen that the DMRG results are not consistent 
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with quasi long-range order. 

In this section we will show that a good qualitative 
description of the DMRG results can in fact be obtained 
from a purely classical treatment of a low-energy effec- 
tive Hamiltonian describing the effects of doping the half- 
filled phase, a phase for which we have seen that the 
basic-Umklapp interaction plays a crucial role. The clas- 
sical model is obtained by neglecting, by hand, the quan- 
tum fluctuations in <f> +p that if kept would have caused 
(j) +p to disorder on length scales ^> 1/6. 



1. Classical model for the SF/DDW phase 

To derive the classical model for the SF/DDW phase, 
we take as our starting point the low-energy effective 
Hamiltonian for this phase at half-filling, which will be 
referred to as !Kq ■ Its form is assumed to be well approx- 
imated by the form for weak interactions, Eqs. (|3.11fl - 
(|3.14() . with values of the effective couplings appropri- 
ate for the system being in the SF/DDW phase for 
intermediate-strength interactions. In order to investi- 
gate the effects of doping, we add to !Ho a chemical po- 
tential term AOi = —f-iN, where 



Af=J2^PXs = - [dxd : 

PXs n 



xV+p 



(5.25) 



measures the number of electrons in the system with re- 
spect to the half-filled case (note that here we work with 
a fixed chemical potential [i instead of at a fixed doping 
S as we have done up until now both in the DMRG and 
bosonization analyses). The fields <f>+ ai Q-p-, and (9_ CT arc 
taken to be locked to values consistent with being in the 
SF/DDW phase, and the resulting gaps will be assumed 
to be large enough that excitations in these fields may be 
safely neglected when discussing the low-energy physics. 

To obtain an effective, classical Hamiltonian from "Hq 
we proceed as follows: (i) Neglect all terms that do not 
contain +p fields, as these can be regarded as constants 
for our purposes, (ii) In the basic-Umklapp term l|3.14f) . 
replace cos2</> +(T , cos2#_ p , cos2</>_ CT , and cos26'_ (T by 
their expectation values (which will generally depend on 
pf). This gives an effective coupling constant for cos 2(f> +p . 
(iii) Neglect the term containing (d x 6+ p ) 2 which gives 
rise to quantum fluctuations in <f+ p . This makes (f> +p a 
purely classical field which allows us to neglect the dis- 
tinction between (f> +p and (4>+p) in the following. Then 
the effective classical Hamiltonian, denoted by 3~C +p , can 
be written 



"X+p — 



dx 



A 



,} +p 



dx 



B cos 2(f> 



+p 



(5.26) 



where A and B are effective positive /^-dependent cou- 
pling constants. Since we are here considering the theory 
on a finite system size L, it should be noted that bound- 
ary effects have not been taken into account in this effec- 
tive Hamiltonian. 



We want to consider a system with a given average dop- 
ing 5 = Nh/2L, where Nh = —N is the number of doped 
holes. Thus AxC is simply a constant —2\/j,\5L which can 
be neglected, and (f>+ p has to satisfy the boundary con- 
dition 



(L) - +p (O) = -ttSL. 



(5.27) 



Since the net spin in the z direction is measured by the 
operator 



2 f L 

= - I dx d x (f +a , (5.28) 



the fact that (d x (f +a ) = in the low-energy subspace 
implies that the number of up- and down-spin electrons 
are equal. Consequently, up- and down-spin electrons 
are removed in pairs by the chemical potential, so that 
Mh is an even integer, and therefore the product SL in 
Eq. (|5.27|) is an integer. Note that this condition is 
satisfied for all figures discussed in Sec. IIVI except Fig. 
Elwhere the presence of an unpaired hole is seen to result 
in disturbances in the density and current patterns. 

The problem we have arrived at consists in finding 
the function cf> +p (x) which minimizes IH_|_„, subject to 
the boundary condition (|5.27|l . That is, the solution 
(f>+ p {x) is the classical ground state configuration of the 
sine-Gordon model in the presence of a finite density of 
solitons)£& equivalcntly, we can think of it as a Frcnkcl- 
Kontorova-like problem^ There is a competition be- 
tween !K+ P which wants to lock <f> +p in one of the minima 
of the cosine potential, and the boundary condition which 
forces <f+ p to have a nonzero slope. 

Once 4> +p (x) has been found, the currents and densities 
are easily obtained from Eqs. (|3.16f) and 15.1|I - H5.2|) . It is 
important to note that since the doping S is determined 
by the term —pAf in the effective Hamiltonian, the 5 = 
form of Eqs. H5.lfH5.2f> must be used. Consequently, the 
staggered rung current can be written 



js{x) = j Q COS(f +p , 



(5.29) 



where jo ~ j-!r(cos(p +a cos0_ p cos6Lo-) ^ 0. The value 
of the amplitude jo will be determined by fitting to the 
DMRG results. Furthermore, the electron density on the 
top and bottom site at rung x are equal and given by 



n(x) 



1 , 1 d(f +p 
7r dx 



(5.30) 



2. Solution of the model 



To solve this problem analytically, we interpret JC+ p as 
the classical action for a particle with mass 2A, located at 
position <f> at time x, with potential energy B cos 2<f> (here 
and in the rest of this paragraph we omit the subscript 
+p on (f)+ p ). It will be convenient to take the potential 
energy to be V{<f>) = B(cos20 - 1) = -25 sin 2 < 
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corresponding to a redefinition of its zero. The sum of 
kinetic and potential energies is conserved, 



(5.31) 



where E is the total energy of the particle. According 
to Eq. H5.27f) . in the time interval L the particle moves 
a distance ttSL, which (since L > 1/5) is bigger than 
7r, the distance between the potential minima. Thus the 
particle is not trapped inside one of the wells of the po- 
tential V(4>), but is unbounded. Hence E > 0, so that 
the velocity d(f>/ dx always has the same sign, which must 
be negative since 4>{L) < 0(0). Eq. I|5.31[) can then be 
integrated to give 



x 



fcF(<t>\-2b), 



(5.32) 



where a = A/E, b = B/E, F((f)\m) is the incomplete 
elliptic integral of the first kind (see Appendix [UJ, and 
xq is an integration constant defined by 4>(xq) = 0. This 
gives 



Xq 



- 2b 



(5.33) 



where am(M|m) is the Jacobian amplitude function. Us- 
ing Eqs. I|5.27[) . H5.33fl . an d (|C4|) . and the assumption 
that L is an integer multiple of 1/5, the parameters a 
and b can be related as 



26K(-2b), 



(5.34) 



where K(m) is the complete elliptic integral of the first 
kind. 

The staggered rung current j s (x) and electron density 
n(x) can then be expressed in terms of the Jacobian el- 
liptic functions cn and dn, respectively: 



n(x) 



Jo cn 



x - xq 



la 



dn 



-2b), 

X - Xq 



2b 



(5.35) 
(5.36) 



Using Eq. (|5.34|) to eliminate a, we see that n{x) can be 
expressed in terms of the unknown parameters b and xq, 
while for j s (x) an additional parameter jo is also needed. 
The parameter xq plays the role of a phase variable which 
distinguishes between different solutions related by trans- 
lation along the x axis, which all have the same energy 
since we have neglected boundary effects. 

Using Eqs. i|5.34[l - i|5.36|) an( i the properties of the ellip- 
tic functions, the following results are easily established: 
n(x) and j s (x) oscillate with wavelengths 1/5 and 2/5, 
respectively, and n(x) is minimal where j s (x) changes 
sign, i.e. at the anti-phase domain walls in the current 
pattern. We emphasize that these properties of the an- 
alytical solution are independent of the actual values of 
the fitting parameters b, j and xq, and are in excellent 
agreement with the numerical DMRG results. 



3. Fits to DMRG results 

In Fig. El we show fits of the analytical expressions 
for j s (x) and n(x) to DMRG results for a 200-rung lad- 
der for 5 = 0.02 (upper panel) and 5 = 0.04 (middle 
panel). The corresponding solutions of (f> +p (x) are also 
shown (lower panel). For both dopings, the model pa- 
rameters, and therefore the DMRG results, are identical 
to those in Fig. ^] The fit is rather good for 5 = 0.04, 
the main deviations occurring near the boundaries due 
to the neglect of boundary effects in the analytical solu- 
tion. However, since by construction both the numerical 
and analytical curves for n(x) integrate to half the to- 
tal number of electrons, deviations between the density 
curves near the boundaries inevitably imply some devi- 
ations also away from the boundaries. These deviations 
come mainly in the form of a horizontal shift between the 
density curves which decreases towards the center of the 
ladder. There is an identical shift between the current 
curves so that the zeros of the current are always at the 
minima of the density. It is as if the boundaries exert a 
"push" on the density and current oscillations, and this 
effect cannot be captured by the analytical solution. 

Similar features are observed in the fit for the lower 
doping 5 = 0.02, but now there are also more pronounced 
differences between the shapes of the analytical and nu- 
merical curves. In the analytical fit the current oscil- 
lations have a more rectangular shape and the density 
minima are deeper. It is natural to attribute these differ- 
ences primarily to the omission of quantum fluctuations 
in the effective model. The step-like nature of the soli- 
tons in the analytical solution for <fi +p (x) becomes more 
pronounced with decreasing doping 5 (see lower panel in 
figure), as the ratio of the soliton separation 1/5 to the 
soliton width becomes larger. We expect that quantum 
fluctuations will be more efficient in smoothing the soli- 
tons the more step-like these solitons are. This picture 
is consistent with the reduced quality of the fit as the 
doping is reduced from 0.04 to 0.02i£& Comparison with 
the DMRG results in Fig. 1171 thus seems to suggest that 
while the solitons in (j) +p (x) would not be much affected 
by the inclusion of quantum fluctuations in the effective 
model for 5 = 0.04, they would be somewhat smoothed 
for <5 = 0.02, leading to less rectangular current oscilla- 
tions and not so deep minima in the electron density. 



4. The CDW and doped D-Mott phase 

The results for the charge density oscillations in the 
CDW and doped D-Mott phase that follow from such a 
classical approach are also in good qualitative agreement 
with the DMRG results. In the CDW phase, the density 
on leg £ of rung x takes the form (obtained from Eq. I|5.2|) 
with (5 = 0) 



n t ,b(x] 



= i+ ld0+p ± c 

TT dx 
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FIG. 17: (Color online) (a) Staggered rung current j a (x) (red 
curves, left y-axis) and electron density per site n(x) (blue 
curves, right y-axis) as a function of rung index x in the 
SF/DDW phase of a 200-rung ladder with 6 = 0.02. Shown 
are DMRG results (diamonds and circles for currents and den- 
sities, respectively, connected by dashed lines to guide the eye) 
and analytical fits (full lines). The model parameters for the 
DMRG calculation are t = V± = 1, U = 0.25, J± = 1.5 and 
Jy = Vji = 0. The values of the fitting parameters used are 
b « 20.0, | jo | ~ 0.61, and x = 0. (b) Plots for the same 
model parameters as above, except that the doping 5 = 0.04. 
The values of the fitting parameters are b w 1.66, |jo| « 0.52, 
and xo — 0. (c) The function <f>+ p (x) corresponding to the 
fits of j s (x) and n(x). 



where n w n 2 fe F (cos cj) +a sin 6L p cos 6_ a ) ^ 0. In the 
doped D-Mott phase the density also has this form ex- 
cept that no = 0. In both phases the solution for 0+ p 
is given by Eq. {533}. as in the SF/DDW phase. Thus 
{AUf, 0) density oscillations with period 1/5, coming from 



the 



jdx term, are predicted to be present in all 



phases (for the CDW phase, the DMRG results show 
that the (2fcp-,7r) density oscillations dominate over the 
(4fcp,0) oscillations in size, however). Furthermore, be- 
cause DMRG suggests that 4>+ p essentially behaves like a 
classical field over the length scales considered here, the 
conjugate field 8 +p is strongly fluctuating, so that DSC 
correlations decay exponentially over these length scales 
even in the doped D-Mott phase. 



VI. SUMMARY AND DISCUSSION 

In this paper we have studied a generalized Hubbard 
model on a two-leg ladder, focusing on a parameter re- 
gion that shows SF/DDW order. The approximate loca- 
tion of this region in parameter space was found from a 
calculation of the phase diagram of the half-filled model 
for weak interactions, which was done by using Abclian 
bosonization and semiclassical considerations to analyze 
the low-energy effective theory resulting from a one-loop 
RG flow. 

The RG/bosonization calculations in this and previous 
worka 37 ! 38 ! 39 suggest that in order to get SF/DDW or- 
der in the generalized Hubbard ladder it is necessary (al- 
though not sufficient) to have the on-site repulsion U be 
less than the nearest neighbor repulsion V± . One class of 
materials that may be promising for achieving a relatively 
small U is organic molecular crystals, due to the large 
size of the highest occupied molecular orbital; bringing 
two of these molecules close together might then conceiv- 
ably give U < V~ Another way of reducing the effec- 
tive U is through an on-site (Holstein) electron-phonon 
interaction^ 

To study (rational) hole dopings 5 away from half- 
filling, finite-system DMRG calculations were done on 
ladders with up to 2 x 200 sites for intermediate-strength 
interactions. It was found that the SF/DDW phase per- 
sists up to quite high dopings (of order 10-20%) and that 
upon doping, currents remain large, with no evidence 
of decay. The rung currents oscillate with wavevector 
(2fcp,7r), corresponding to oscillations with wavelength 
2/5 in the staggered rung current. The currents coexist 
with (4fcp, 0) charge density oscillations, with two doped 
holes per wavelength 1/(5, which also are not found to 
decay. 

The factor of two ratio between the periodicities of 
the staggered rung current and the rung charge density, 
and the related fact that the maxima of the hole density 
are located at the anti-phase domain walls of the rung 
current, imply that the SF/DDW phase is an example 
of a phase with "topological doping."— In this respect, 
the rung current in the SF/DDW phase is the analogue 
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of the spin density in the so-called "stripe phases" in 
doped antiferromagnets, which have coexistence of spin 
and charge density order ^ 

Charge density oscillations which do not seem to de- 
cay are also found in the neighboring CDW and doped 
D-Mott phases, in which currents decay exponentially. 
Both these phases have (4fci?,0) charge density oscilla- 
tions; the CDW phase in addition has (2kp,ir) charge 
density oscillations. DSC correlations decay exponen- 
tially in all phases (except possibly deep inside the doped 
D-Mott phase where the numerical results for the decay 
are more difficult to interpret). 

We have shown that most of these DMRG results for 
the doped ladder can be qualitatively understood from 
weak-coupling RG/bosonization arguments. However, 
there is one apparent exception: The DMRG results are 
consistent with true long-range order in the currents in 
the SF/DDW phase and true long-range order charge 
density order in all phases. As far as we know, such a 
true long-range order scenario is not forbidden by any 
exact theorems that are directly applicable to the lat- 
tice Hamiltonian studied here; in particular, the Mermin- 
Wagner theorem is respected. On the other hand, the 
RG/bosonization arguments imply that such a scenario 
would require the symmetric charge field c/)+ p to be- 
come locked and that this can only happen if multiple- 
Umklapp interactions are relevant in the RG sense. How- 
ever, a calculation of the scaling dimension of these inter- 
actions suggests that they are irrelevant for the dopings 
considered, and thus the true long-range order scenario 
does not appear to be supported by the RG/bosonization 
arguments. These arguments further suggest a crossover 
from non-decaying current/charge density correlations to 
power-law decay at a length scale of order 1 /8, due to the 
basic-Umklapp interaction being unimportant on much 
longer length scales^ For the same reasons, in the doped 
D-Mott phase a crossover from exponential to power-law 
decay in the DSC correlations would be expected at the 
same characteristic length scale. 

In light of this discrepancy between the predictions of 
the DMRG and RG/bosonization method for the long- 
range order question, a detailed scrutiny of both methods 
might be necessary, to see if one could possibly identify 
an omission or shortcoming of the conventional analysis 
that, if corrected, would resolve the discrepancy. While 
any investigations along these lines are beyond the scope 
of the present work, we will in the next paragraphs com- 
ment on the issues that might appear most salient. 

For the DMRG method, a standard concern is whether 
boundary and/or finite-size effects could cause the sys- 
tem to appear to have true long-range order even if there 
in reality is none. However, the boundaries are in fact 
found to suppress SF/DDW order, not enhance it. Fur- 
thermore, the result of a finite-size test for Luttinger liq- 
uid behavior (quasi long-range order), that involved the 
evolution of the mid-ladder density with ladder length, 
was negative. Thus there seems to be no reason to ex- 
pect the DMRG results to change qualitatively for larger 



systems than we have studied here. The DMRG calcula- 
tions are also "internally consistent" in the sense that the 
DSC correlations in the doped D-Mott phase are found 
to be consistent (or at least not inconsistent) with expo- 
nential decay, in agreement with the bosonization result 
when <p +p is locked. 

As the true long-range order scenario implies that <fi +p 
is gapped, it would have been desirable to use DMRG to 
calculate the charge gap directly, as function of ladder 
length, and see if it extrapolates to a nonzero value in 
the thermodynamic limits However, this calculation is 
very difficult and we have unfortunately not been able to 
do it. 

It has been known for a long tim o 71 i 92 that DMRG 
works with ansatz states of the form of matrix prod- 
uct statesi^&Si Within this class of ansatz states it finds, 
with some qualification of no importance to the present 
argument, the variationally optimal state as approxima- 
tion to the true ground state. One might think that this 
methodological bias might sometimes algorithmically fa- 
vor a state that is not a good approximation to the true 
ground state (in the present case, suggesting long-range 
instead of quasi long-range order). To our knowledge 
there is however no known case where converged DMRG 
calculations have seriously misrepresented the true quan- 
tum state where it is known from other methods. More- 
over, matrix product states would typically bias more in 
favor of shorter-ranged correlations. Yet, this scenario 
cannot be totally excluded. 

One concern about the RG/bosonization method is to 
what extent its predictions can be trusted when the in- 
teractions are not weak. Wu et al. have recently argued 
that the quasi long-range order predictions would be ex- 
pected to be valid also for the interaction strengths used 
in the DMRG calculations presented here although the 
parameters in the low-energy effective theory might be 
strongly renormalizedi2£ A concrete example of the va- 
lidity and usefulness of RG/bosonization arguments for 
analy zing a ladder with strong interactions was given in 
Ref. |45l In particular, for 3/8 filling in the t-J lad- 
der with J/t ss 0.25, DMRG calculations showed ev- 
idence for a charge density wave with true long-range 
order; the values for K +p extracted from the DMRG 
results for J/t > 0.25 were consistent with the critical 
value K+ p = 0.125 predicted by RG/bosonization for 
the onset of true long-range order4^ But this example 
does of course not imply correctness of the weak-coupling 
RG /bosonization predictions in the model considered by 
us. 

Another, much more speculative concern is whether 
some of the RG arguments used in this paper might break 
down in a more fundamental way. A scenario for the 
breakdown of the perturbative RG at certain quantum 
critical points has recently been discussed in Ref. |95l The 
breakdown occurs because the expansion coefficients of 
terms that appear in the RG equations at two-loop order 
depend singularly on a dangerously irrelevant variable, 
which causes them to diverge. Physically, the effect is 
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related to the presence of another diverging time scale in 
addition to the critical one. For coupled one-dimensional 
chains, nonlinear corrections to the dispersion (i.e., band 
curvature terms), which like the Umklapp interactions 
reflect the underlying lattice structure of the problem, 
do seem to behave in a way that might be characterized 
as dangerously irrelevant with respect to some transport 
properties (Coulomb drag)i2Si2i Whether as a result of 
this something akin to the breakdown described in Ref. 

might occur in the model studied here is an intriguing, 
but at this stage highly speculative, question. 

Although resolving the long-range order question is an 
important theoretical issue, as a practical matter it may 
make little difference whether there is true long-range 
order or only algebraic order that decays very slowly. 
Physical realizations of systems such as we consider in 
this paper will either involve weakly coupled ladders, or 
a full two-dimensional lattice. If isolated ladders can only 
show quasi-long-range ordered current and charge density 
correlations, even extremely weak coupling between such 
ladders could suffice to stabilize true long-range order. 
We note that recent results of inelastic neutron scattering 
measurements on Sri4Cu2404i have been interpreted as 
a possible signature of orbital currents^ Clearly it would 
be very interesting if further experiments on this material 
were to corroborate this interpretation. 



order parameters for the SF/DDW and CDW phases, 
Eqs. (|3.15|) - (|3.17|) . contain the two anti-Hermitian op- 
erators h s = r)2 S r]is, s =T,j, while the bosonized order 
parameters for DSC and SSC correlations, Eqs. (|3.19|) - 
(|3.20() . contain the two anti-Hermitian operators h' x = 
7y AT ^ Ai; A = 1,2. From T 2 = -h 2 s = -hf = 1, it follows 
that r has eigenvalues ±1, while {h s } and {h' x } have 
eigenvalues ±i. 

Since [T,h s ] = [T,h' x ] = [h s ,h s >] = [h' x ,h' x ,] = 0, it is 
possible to simultaneously diagonalize the order parame- 
ter of a given phase and the Hamiltonian in the Klein fac- 
tor space. (Note that although [h s ,h' x ] ^ 0, this causes 
no problem for such a simultaneous diagonalization, be- 
cause each order parameter contains either {h s } or {h' x }, 
not both.) 

A consistent choice of eigenvalues can then be deduced 
from the relations T = —h^h^ = h\h!^. In this paper we 
choose the eigenvalues T = 1, h a = i, and h' x = i(—l) x . 



APPENDIX B: RENORMALIZATION GROUP 
EQUATIONS AND INITIAL CONDITIONS FOR 
CONTINUUM COUPLINGS AT HALF-FILLING 

The one-loop RG equations can be derived using the 
operator product expansion (OPE)i£2iffii§2i2£ At half- 
filling they read 
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APPENDIX A: KLEIN FACTOR CONVENTIONS 

Here we explain the conventions used for the Majorana 
Klein factors^i^ in the bosonized versions of the Hamil- 
tonian and the various order parameters considered in 
this paper. 

The nonquadratic part of the Hamiltonian density, 
H ( j b) +nf ] in Eqs. (|3~T^|) - (|3~TIjl . contains the Hermitian 
operator T = ^lfTyix^T^i- Furthermore, the bosonized 
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Sip = fJt P + Yg-9t<T — Stui — 9tu2 ~ 9tui9tu2, 

2^ 2 2 2 2 

9xp = ~9tp — ~l6^ ta ~ ^ XU ~ ^ tUl " ^ tu2 ~ 9tul9tu2, 



9 In 



-9ia - -j9t<r + 2 9t P 9ta ~ 4g tul - 4g tul g tu2 , 



9xo- = ~9xo- ~ - 2 9t P 9to- ~ 4.g tu2 - ^9tu\9tu2, 

3 

9t P = 2gt P {gip - gxp) + -gto(gio - sw) 

— 9xu(gtul — 9tu2), 

gto- = 2gt P (gio- - g X o) + gto(2gip - 2g xp - gio - gxo-) 
+ ^gxu{gtui + gtu2), 



9xu — ~^9xp9xu — gtui ( 2gt P — -^gto 



+ gtu2 ( 2gt P + 2 gta J ' 

3 1 

gtui = —gtui [ 2gi p + 2g xp + -^9ia — -^9xo ) — giogtu2 

1 



+ gxu ( -2gt P + -^gtoj , 

1 3 

gtu2 — —gtu2 [ 2<7i p + 2g xp — 2# lcr + i^9xo ) — gxogtui 



cm ( 2gt P + -^9to 



(Bl) 



Here iji = 2itv Fdgi/dl, where / is related to the running 
cutoff scale as a(l) = ae l . These RG equations were first 
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derived in Ref. [5£ 
are found to be 



The initial values for the couplings 
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2 2 



-V[i - -Ji - -Jii 



-c/ - -vi — VTi — - Ji 

2 2 11 8 



-Ji 



(B2) 



Eqs. (fBlT) and (|B2|) hold for t x = t, i.e. fc F2 = 2&fi = 
27r/3. As a partial check of the internal consistency of 
these equations, it is instructive to generalize them to 
arbitrary values of t±/t (though < 2 so that both bands 
are occupied) and by also allowing for density-density 
and spin exchange interactions along the plaquette di- 
agonals with strength V and J', respectivelyiiffi It can 
then be shown that the generalized equations have the 
correct symmetries in special limiting cases: U(4) sym- 
metry when t± = 0, U = V±, V\\ = V , J\\ = J± = J' = 0; 
U(2)xU(2) symmetry when t± = Vj_ = V = J± = J' = 
(i.e. independent legs), and SO (5) symmetry whe n 58 i 68 
Ji = 4(C7 + Vj_), Vj| = V = J,, = J' = 0. (Note that 
when tj_ = three new continuum couplings are allowed 
in addition to the nine present for a generic t±.) 



APPENDIX C: ELLIPTIC INTEGRALS OF THE 
FIRST KIND AND ASSOCIATED JACOBIAN 
ELLIPTIC FUNCTIONS 

In this Appendix we summarize notationifli for and 
some basic properties of the elliptic integrals of the first 
kind and the associated Jacobian elliptic functions that 
are encountered in the analytical solution of the Prenkel- 
Kontorova-like problem in Sec. IVE2I 

The definitions of the incomplete and complete elliptic 
integrals of the first kind are, respectively, 



F{(t>\m) = 



4> 



da 



'o \J 1 — m sin 2 a 
K(m) ee F(tt/2|to). 



(CI) 
(C2) 



We will only be concerned with negative values of m here. 
The inverse function of u = F(d>\m) exists and is defined 



4> ee am(u|m). (C3) 
This function is known as the amplitude for the Jacobian 
elliptic functions. It is odd in it, and satisfies 



am(u + 2 A' (to) | to) = &m(u\m) + 7T. 



(C4) 



The Jacobian elliptic functions encountered in our prob- 
lem are 



cn(w|m) 



COS( 



(C5) 



dn(u\m) — — - = \/ 1 — to sin 2 <p. (C6) 
du 



These two functions are even and periodic in u with pe- 
riod 4 A (to) for cn and 2 A (to) for dn. 
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